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ABSTRACT 

We report the probable identification of the X-ray counterpart to the 7-ray 
pulsar PSR J2021+4026 using imaging with the Chandra X-ray Observatory 
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ACIS and timing analysis with the Fermi sateUite. Given the statistical and sys- 
tematic errors, the positions determined by both satellites are coincident. The 
X-ray source position is R.A. 20^21"30".733, Decl. +40°26'46.04" (J2000) with 
an estimated uncertainty of l."3 combined statistical and systematic error. More- 
over, both the X-ray to 7-ray and the X-ray to optical flux ratios are sensible 
assuming a neutron star origin for the X-ray flux. The X-ray source has no 
cataloged infrared-to-visible counterpart and, through new observations, we set 
upper limits to its optical emission of i' > 23.0 mag and r' > 25.2 mag. The 
source exhibits an X-ray spectrum with most likely both a powerlaw and a ther- 
mal component. We also report on the X-ray and visible hght properties of the 
43 other sources detected in our Chandra observation. 

Subject headings: pulsars: individual (PSR J202H-4026) — supernova remnants: 
individual (SNR 78.2+2.1 = 7-Cygni SNR) — X rays: general 

1. Introduction 

Pulsars are often depicted as relatively simple objects: a highly magnetized, fast rotating 
neutron star (NS) emitting radiation along its poles. Most emission models start from this 
basic picture, in principle common to all isolated pulsars. The rotation of the magnetic 
dipole induces an electric fleld near the polar caps in vacuum, which is then short circuited 
by an e^ pair plasma (Ruderman & Sutherland 1975), or in space-charge limited flow, which 
is screened at a pair formation front (Arons & Scharlemann 1979). Models seek then suitable 
regions of particle acceleration. Although the radio emission is thought to originate in the 
polar regions of the magnetosphere (e.g Michel 1991 and references therein), the high-energy 
emission is now thought to originate from the outer magnetosphere, with the original polar 
cap concept (Daugherty & Harding 1996) having been superseded by outer-gap (Cheng, 
Ho & Ruderman 1986; Romani 1996) and slot-gap (Arons 1983, Muslimov & Harding 2004) 
models. The picture that has emerged is that radiation is emitted in a continuous, very broad, 
spectral range, with curvature radiation producing most of the GeV emission (Romani 1996) 
and synchrotron and/or Compton scattering by cascade products producing the near-infrared 
to soft 7-ray emission (Takata & Chang 2007, Harding et al. 2008). 

The high energy emission properties of pulsars were revealed in the 1990s when use of 
the EGRET experiment on the Compton Gamma-Ray Observatory (e.g. Thompson, 2001 
and references therein) led to the multiwavelength spectral characterization of half a dozen 
7-ray pulsars, including the discovery of the radio-quiet pulsar Geminga (Halpern & Holt 
1992), the second brightest steady GeV source in the 7-ray sky. Amongst the EGRET 
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legacy was a sample of 170 unidentified sources, 74 of which are at |6| < 10° (Hartmann et 
al. 1999). It has now been found that ^ 43 EGRET unidentified sources have counterparts 
in the Fermi Large Area Telescope (LAT) first year catalog (Abdo ct al. 2010b; IFGL) and 
a large fraction of those lying on the Galactic plane turned out to be pulsars (Abdo et al. 
2008, 2009a, 2009b, 2010a, 2010b), a result anticipated by several authors (Yadigaroglu & 
Romani 1995, Cheng & Zhang 1998, Harding & Muslimov 2005). 

In fact, it is very plausible that many more radio-loud or radio-quiet pulsars are hidden 
in the unidentified Galactic LAT sources, although estimates are highly dependent on details 
of the different emission mechanisms. At this point it is still unclear what makes some 7- 
ray pulsars radio-loud and some radio-quiet. Understanding the distinct properties of the 
individual sources will surely lead to a better understanding of models for the emission 
mechanisms, for example, the connection between the radio and the near-infrared to the 
7-ray component. 

The source 2CG 078-1-2 was one of about twenty 7-ray sources detected by the COS-B 
satellite 30 years ago (Swanenburg et al. 1981). The source is in the fine of sight towards the 
supernova remnant G78. 2-1-2.1. The remnant comprises a 1-degree-diameter, circular, radio 
shell with two bright and broad opposing arcs on its rim (Higgs, Landecker & Roger 1977, 

Wcndkcr, Higgs & Landecker 1991). The remnant has a kinetic distance of about 1.5 kpc 
(Landecker, Roger, and Higgs 1980; Green 2009) and is estimated to have an age of 5400 yr 
(Sturner & Dermcr 1995). 2CG 078+2 is often cited as 7-Cygni due to its proximity to the 
second magnitude foreground star (niy = 2.2, spectral type F8Iab) that lies on the eastern 
edge of the remnant although there is no other physical association between the supernova 
remnant (SNR) and the star. A small HII region, located close to the star, forms the 7-Cygni 
nebula. 

The 7-ray source, next cataloged as 3EG J2020+4017 (Hartmann et al. 1999), was 
suspected to be either extended emission from the SNR or a pulsar formed in the supernova 
event (Sturner & Dermer 1996). Soon after its June 2008 launch, the Fermi 7-ray Space 
Telescope highhghted the discovery of a radio-quiet pulsar in the CTA-1 SNR (Abdo et 
al. 2008) as a first light result. This was followed by the discovery of twenty-three other 
7-ray pulsars using blind search techniques (Abdo et al. 2009b, Saz Parkinson ct al. 2010), 
among them PSR J2021+4026 lying within the EGRET error circle of 3EG J2020+4017. 
PSR J202 1+4026 is a 265-ms pulsar with a spin-down age of 77 kyr and a total spin-down 
luminosity of 1.1x10^^ erg/s. Interestingly, most of the pulsars discovered by the LAT 
with blind searches have not been seen at other wavelengths. In fact, only three of the 26 
discovered to date have been detected in radio (Camilo et al. 2009; Abdo et al. 2010c) 

Even before the Fermi discovery of pulsations, and because the EGRET 7-ray proper- 
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ties of the source matched those of known 7-ray-pulsars (hard, steady and in the Galactic 
plane), the search for the radio and X-ray counterpart of this source began (Brazier et al. 
1996; Becker et al. 2004, Weisskopf et al. 2006). No radio pulsar is associated with PSR 
J2021+4026 with searches having now been conducted down to lOOyuJy at 1665 MHz (Trepl 
et al. 2010), 40 /xJy at 820 MHz (Becker et al. 2004), and to 11/iJy at 2 GHz (Ray et al. 
2011). There is an extended radio source (GB6 J2021+4026) in the vicinity which appears 
to be positioned more or less symmetrically (see Figure 8 of Trepl et al. 2010) about our 
best position for the X-ray counterpart but there is no evidence that it is associated with 
the 7-ray pulsar. 

We have previously (Becker et al. 2004, Weisskopf et al. 2006) performed Chandra 
observations (ObsIDs 3856 & 5533) aimed at different portions of the 7-Cygni field and these 
pointings were based on the best known 7-ray positions available at the time. The latter 
Ciiaudra observation, ObsID 5533, overlapped with about 3/4 of the current Fermi-LAT 99%- 
confidence positional error circle and detected several potential X-ray counterparts, including 
the source designated as "S21" as reported in Weisskopf et al. (2006). Subsequently, Trepl 
et al. (2010) reanalyzed the available Chandra data and also used XMM-Newton data to 
search for a counterpart. They identify source 2XMM J202131-I-402645, at virtually the 
same location as the Chandra source "S21", as the likely counterpart because the X-ray 
source position fell within the ~4'-radius Fermi-LAT 95%-confidence positional error circle 
at the time (OFGL, Abdo et al. 2009b). We note that the most recent Fermi-LAT error circle 
radius (2FGL, Abdo et al. 2011) is 0.6' (§3) and no longer includes "S21". We re-evaluate the 
situation using both the greatly improved LAT source localization and the position measured 
from timing of PSR J2021+4026. 

The current work describes a new observation that completes the Chandra coverage of 
the field. As a result of these observations, we conclude that the source originally labelled 
"S21" by Weisskopf et al. (2006), 2XMM J202131+402645 by Trepl et al. (2010), and 
source #20 in the new observation, remains the most probable X-ray counterpart to PSR 
J2021-(-4026. Moreover, we also show that the X-ray source is dominated by thermal, not 
powerlaw, emission. In this regard, it is interesting to compare this source with Geminga, 
one of the best studied radio-quiet 7-ray pulsars, which is older (340 ky) and less luminous 
(~ 3 X 10^^ erg/s) than PSR J2021+4026 but has a similar period and period derivative. 

§2.11 describes the analysis of the X-ray image; §2.21 describes the analysis of the X-ray 
spectrum, especially of the X-ray source #20; §2.31 compares the X-ray spectral properties 
of source #20 to Geminga and CTA 1; §3] describes new Fermi-LAT pulsar timing of PSR 
J202H-4026; §1] describes our search for an optical counterpart; and ^discusses properties 
of the other Chandra sources in the field. We provide summary conclusions in ^ 
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2. Chandra X-Ray Observations and Data Analysis 

We obtained a 56-ks Chandra observation (ObsID 11235, 2010 August 27) using ttie 
Advanced CCD Imaging Spectrometer (ACIS). Here we report the data taken witli the 
back-illuminated CCD ACIS S3 in faint, timed-exposure mode with 3.141-s frame time. 
Background levels were nominal throughout the observation. Standard Chandra X-ray Cen- 
ter (CXC) processing provided accurate aspect determination. 

Starting with level-1 event lists, we reprocessed the data using the CIAO v4.2 tool 
acis_process_events to remove pixel randomization in order to improve the on-axis point 
spread function (PSF), thus enhancing the source-detection efficiency and positional accu- 
racy. In searching for sources, we utilized events in pulse-invariant channels 35-550 corre- 
sponding to 0.5 to 8.0 keV. 

To verify the Chandra position accuracy, we compared the pointing parameters (RA, 
Dec, roll) given by the Chandra data-processing to that using the 19 2MASS sources (Sec- 
tion |53l2]) with high probability X-ray counterparts. For these calculations we used 0.075" 
per axis for the 2MASS position uncertaint}0. Assuming no change to the Chandra param- 
eters yielded an acceptable positional fit for these 19 sources with of 24 for 38 degrees 
of freedom. Allowing the three pointing parameters to vary marginally improved to 21 
for 35 degrees of freedom and would imply the following corrections: (0.14 ± 0.10)" (RA), 
(—0.04 ± 0.08)" (Dec), and (43 ± 78)" (roll). However, these corrections, even allowing for 
errors, being negligible we did not correct the Chandra positions. 

2.1. X-Ray Image Analysis 

We searched for point-like X-ray sources employing techniques described in Tennant 
(2006), using a circular- Gaussian approximation to the PSF and setting the signal-to- noise 
(S/N) threshold for detection to 2.4. The resulting background-subtracted point-source 
detection limit is about 7 counts, with fewer than 1 accidental detection expected over the 
field. Based upon tests on Chandra deep fields, this approach finds all X-ray sources in those 
fields down to 10 counts, which we thus regard as the completeness limit. Figure [T] shows 
the ACIS-S3 image with a small circle at the position of each Cliaiidra-detected source. 

Table [1] tabulates the X-ray properties of the 44 Ciiandra-detected sources, with the 
source number in column 1. Columns 2-5 give, respectively, right ascension RA, declination 



http://www.ipac.caltech.edu/2niass/releases/allsky/doc/explsup.htnil 
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Dec, extraction radius 6'ext, and approximate number of X-ray counts Cx detected from the 
source. Column 6 lists the single-axis RMS error = [(cTp^p/Cx) + c^sys]"^''^ ^^e X-ray- 
source position, where cxpsF is the dispersion of the circular Gaussian that approximately 
matches the PSF at the source location and CTgys is a systematic error. Uncertainties in the 
plate scal^ imply cTsys ~ 0'.'13: To be conservative, we set cxsys = 0'.'2 (per axis). Column 7 
gives the radial uncertainty egg = 3.03 ax in the X-ray position — i.e., xi = 9-21 = 3.03^ 
corresponds to 99% confidence on 2 degrees of freedom — for inclusion of the true source 
position. Columns 8 & 9 are color ratios defined and discussed in § 15. 1[ 

In view of the spin-down age and energetics of PSR J2021+4026, the possibility exists 
that a pulsar wind nebula (PWN) may also be present in the X-ray image of the field. 
We searched for moderately extended sources in the field and identified three features of 
interest. Two are located near the edge of the S3 CCD which increases the likelihood of 
their being false positives. The third feature lies ~6" west of source #20. Simultaneously 
fitting the combination of this feature and source #20 to circular Gaussians (plus a constant 
background) results in a Gaussian width of 2.2" and a total of 22 X-ray events within this 
extended feature. 

We also compared the spatial distribution of events in and around source #20 to a model 
of the Chandra/ ACIS PSF valid for its location relative to the aimpoint and its characteristic 
energy using the PSF library psflib v4.1 (CALDB v4.2). Source #20 is consistent with being 
point like. 

2.2. X-Ray Spectral Analysis 

We used the XSPEC (v. 12. 5) spectral-fitting package (Arnaud, 1996) to perform spectral 
fits to the 44 X-ray point sources. We treat source #20 separately in §2.2. 1[ Data were 
binned to obtain at least 10 counts per spectral bin before background subtraction. The 
background was determined after masking off a circular region around each of the 44 detected 
X-ray sources corresponding to a circle of radius 25 times the uncertainty listed in column 
6 of Table [1] and then averaging over the remaining pixels. Individual response matrix and 
ancillary response files appropriate to each source position were created for this analysis 
using the mkacisrmf and mkarf tools available in ClAO version 4.2. 



http://cxc.harvard.edu/mta/ASPECT/aca_plate_scale/ 
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2.2.1. The X-Ray Spectrum of Source #20 

A sufficient number of counts were detected from source ^20 to perform more exten- 
sive spectral analysis. Figure [2] shows the X-ray spectrum of source #20 with background 
subtracted. Data were again binned to obtain at least 10 counts per bin before background 
subtraction. The background is less than 3% of the signal plus background from the region 
that includes source #20. 

We begin our spectral analysis by ffist considering single- component spectral models 
(Table [2]) with a multiplicative absorption component: an absorbed power law (powerlaw 
in XSPEC); an absorbed black body (bbodyrad); and three different absorbed neutron star 
atmosphere modelsjfl namely, nsa (Pavlov et al. 1995), nsmax-1260 and nsmax-130190 (Ho, 
Potekhin & Chabrier 2008). 

For computing absorption, we utilized abundances (XSPEC's wilm) from Wilms, Allen, 
& McCray (2000) with cross-sections (vern) from Verner et al.(1993) and allowed for inter- 
stellar extinction by grains using the model (tbabs) of Wilms, Allen, & McCray (2000). 

All of these models provide statistically adequate ffis to the data in the 0.5-8.0 keV 
range (Table [2]). (We note that so did a ffi to an absorbed mekal model.) From the Hi in 
the Galaxy (Dickey & Lockman, 1990), one infers a column density A^h ~ 1-4 x 10^^ cm~^ 
through the Galaxy in the direction of source #20, implying that values below this should be 
expected if there is no circumstellar absorption, and this is indeed the case (Table [2]). Next 
we posit that the very steep power law index of almost 5 is not physical, but indicative of a 
soft, thermal component. The bbodyrad model's normalization, (-RkmZ-^iokpc) ~ 0.80, where 
-Rkm is the radius of the emitting area in km and -Diokpc is the distance to the source in units 
of 10 kpc, corresponds to an emitting area of i?km = 0.13, assuming the association of the 
pulsar with G78. 2+2.1 at a nominal distance of -Diokpc = 0.15. This is much smaller than a 
typical NS radius, but not very different than the standard polar cap radius, ^ RjsissinO ~ 
R%{2ti/cPY/'^ which for P = 0.265 s and Rns = 13 km is 0.42 km. 

Conversely, if we assume the emission comes from a NS with a characteristic radius of 
13 km, then the bbodyrad norm implies a source distance of ~ 145 kpc, well outside the 
Galaxy. This distance estimate drops for the different NS-atmosphere models, e.g. down 
to 13 kpc (Table 12]), still somewhat distant to remain within the Galaxy. In addition, the 
temperature estimates of the (cooling) NS, range from log(Too) > 6.5 to 6.1, with the precise 
value depending on the model (Table [2]). These estimates are consistent with those expected 



■^A NS radius of 12.996 km and mass of 1.358 Mq is assumed throughout for purposes of computing 
effects of gravitational redshift in the neutron star atmosphere models 
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for a pulsar of an estimated age somewhere between 5400 and 77000 years, depending of 
course on the equation of state, composition of the heat-blanketing envelope, and the degree 
of superfluidity in the star's core. (See Yakovlev et al. 2010 and references therein for recent 
details on the subject of cooling NSs.) 

We then ask whether or not combining a power law with the other models is indicated 
by the data. Table [3] tabulates the change in x^, the f-statistic, and the probability that 
combining models has significantly improved the quality of the fit. The table also tabulates 
he derived spectral parameters. In all cases combining various thermal models with a power 
law does improve the quality of the fit with a confidence better than 2-sigma, but not 3-sigma. 
Moreover, these two-component models allow a wide latitude for the uncertainties of the best- 
fit parameters and hence these parameters are not as constrained as one might wish. This 
follows from the fact that all of the single models (Table themselves provide statistically 
adequate fits. Thus, the 3-sigma contours for any two-component model allows for one or 
the other of the component models to have a zero norm - e.g. a power law component is not 
completely required by the data at this level of significance. Keeping this proviso in mind, 
we continue to examine the two-component models. Certainly the physical interpretation 
of the data is perhaps more sensible when both components are introduced, especially in 
light of the strong 7-ray emission which cannot arise from any single-component "thermal" 
model spectrum. In addition, both the inferred NS temperatures and distances implied by 
the two-component models (Table [3]), with the possible exception of the bbodyrad+powerlaw 
model, are consistent with what one might expect for a young cooling neutron star within 
the galaxy. 



2.2.2. Discussion of the Spectrum of Source #20 

The considerations of the single-component models in the previous section lead one to 
conclude that we are perhaps seeing a hot spot of size comparable to the polar cap rather than 
thermal emission from the entire surface of the neutron star and with a temperature higher 
than one expects from cooling of a neutron star at an age > 5000 years. It is thus possible 
that the bulk of the emission comes from heating of the polar cap by backflowing accelerated 
particles. The expected luminosity and temperature from heating by positrons produced by 
curvature radiation of primaries in a space-charge limited flow model is ~ 10^^ ergs~^ and 
log(T+) ~ 6.3 (Harding & Muslimov 2001). This luminosity is deposited over an area roughly 
that of the polar cap and radiated on a timescale less then the heat diffusion timescale across 
field lines to other areas of the NS. The temperature T+ is close to that determined from 
the bbodyrad and nsa model fits allowing for the gravitational redshift. It is therefore quite 
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possible that the X-ray emission from PSR J2021+4026 is dominated by polar cap heating. 

On the other hand, many well-studied neutron stars that exhibit both X-ray and 7- 
ray emission have composite X-ray spectra, showing non-thermal power law magnetospheric 
emission and/or hot polar cap emission in addition to the lower-temperature full surface 
thermal emission (e.g. Geminga, PSR B0665+14, PSR B1055-52, see De Luca et al., 2005). 
While the present data quality do not demand such two-component models, the inferred 
distances from such two-component fits, e.g., ~ 6 kpc for the NS atmosphere models (Table 
[3]) become reasonable for plausible stellar radii assuming full surface emission. This distance 
is larger than the kinematic distance to the SNR 78.2-1-2.1 of 1.5 kpc but is comparable to 
the distance to the Cygnus arm at Galactic longitude ~78°. Thus, it is possible that a two- 
component model, with full-surface cooling and magnetospheric power law emission present 
at lower levels, might be needed to accurately describe the emission physics. In practice, 
decomposing such complex X-ray spectra requires good statistics and phase-resolved X-ray 
spectroscopy. For a target this faint, extremely long observations or next-generation X-ray 
satellites are clearly required. 



2.3. Comparison of the Spectrum of #20 to Geminga and CTA 1 

There are similarities and differences between the X-ray spectrum of PSR J2021-I-4026 
and two of the other radio-quiet 7-ray pulsars with detected X-ray emission, Geminga and 
CTA 1. First PSR J2021+4026, like Geminga (Jackson and Halpern, 2005) and CTA 1 
(Caraveo et al. 2010), may also be characterized by two spectral components, a thermal 
component and power law. In Geminga, however, the power law component begins to 
dominate above about 0.5 keV and log(Too) is about 5.7. For PSR J2021-I-4026 the power 
law dominates above 2.5 keV and log(Too) is higher, as one would expect as PSR J2021-I-4026, 
based on its spin-down age, is younger. 

The spectral indices for the power law components for source #20 and Geminga are not 
dissimilar, but one needs to recognize the large uncertainty in the measurements reported 
here. Another, possibly important, spectral difference between the two X-ray spectra is that, 
in the case of Geminga, the blackbody component gives an emission radius that is plausible 
for a NS radius. This is not so for PSR J2021-I-4026. In this case, the neutron star atmosphere 
models seem to yield more physically reasonable parameters than the bbodyrad+powerlaw 
model. If we assume that the younger and hotter star still has an atmosphere while the older 
Geminga does not, then these results are sensible. Finally, there is a weak extended emission 
feature near source #20 that may be indicative of a PWN. If so, then it extends no more 
than 0.04 to 0.17 pc from source #20 (assuming a distance of 1.5 to 6.0 kpc, respectively) 
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and contributes ~7% of the X-ray counts detected from source #20 and its surroundings. 
Emission associated with the Geminga PWN has a similar extent and contributes 10% of 
the non-thermal X-ray flux of the pulsar, but only about 1% of the total flux (Pavlov et al. 
2010). 

For CTA 1, the measured temperature, powerlaw index and emission radius are log(Too) = 
6.08, r = 1.3 and r = 0.64 km for a powerlaw+bbodyrad model, and log(Too) = 5.78, 
r = 1.25 and r = 4.92 km for powerlaw+nsa, with slightly lower ior the former (Caraveo 
et al. 2010). In both cases, the emitting radius is significantly smaller than a standard NS 
radius. Thus, similar to PSR J2021+4026, CTA-1 shows a possible heated polar cap com- 
ponent, with a temperature very close to the model prediction of log(Too) — 6.2 (Harding 
& Muslimov 2001), however, the emitting radius in this case is a factor of 2.5 larger than 
the polar cap radius. CTA 1 does not show evidence for a cool component, with the upper 
limits making it unusually cool for its age. 



3. Fermi-LAT Localization and Timing Analysis 

The Fermi-LAT normally localizes 7-ray sources using the incident 7-ray photon di- 
rections. The LAT has a point spread function that is strongly energy dependent, with a 
resolution of about 0.8° at 1 GeV. For a bright source, however, locahzation with arcminute 
accuracy is possible. The source in the second LAT catalog (Abdo et al. 2011) that corre- 
sponds to PSR J2021-F4026 is 2FGL J2021.5+4026. The catalog position for this source is 
R.A. 20'*21'"34M, Decl. +40°26'28" (J2000), with a 95% confidence radius of 0.60'. 

For pulsars, one can use timing techniques better to localize the source, independent of 
the photon direction measurements, as described in Ray et al. (2011). The position deter- 
mination from timing of this pulsar is hampered by the large contribution from rotational 
instabilities common in young pulsars and manifest as "timing noise". Therefore, wc have 
taken two different approaches to try to confirm the association between PSR J2021+4026 
and source ^^20. For this analysis, we used 7 rays detected by the LAT from 2008 Aug 4 
to 2011 Mar 12, selecting only those within 0.8° from the previous best position (Ray et 
al. 2011), and with energies greater than 400 MeV. These "cuts" were chosen to maximize 
the significance of the pulsation. We chose only photons belonging to the most restrictive 
"diffuse" class according to the "Pass 6" instrument response functions (see Atwood et al. 
2009), which have the lowest background contamination. Furthermore, we selected only pho- 
tons with a zenith angle of < 105°, to reduce contamination due to secondary- atmospheric 
7 rays. 
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As a first test, we evaluate the significance of the 7-ray pulsations by assuming the 
pulsar is at each of the candidate X-ray source locations seen in our Chandra observation 
(Table [T]). For each candidate location, we transform arrival times to the barycenter using 
the X-ray position and then use the prepfold routine from the PRESTO pulsar package 
(Ransom, Eikenberry & Middleditch 2002) to find the frequency (/), and its first and second 
derivatives, / and /, which maximize the statistical significance of the pulsation. Figure [3] 
shows the results from this exploratory search, where it is clear that source #20 gives the 
maximum significance for pulsation using this algorithm. This indicates that of the possible 
X-ray sources in the field, source #20 is the most likely X-ray counterpart. 

Next, we use pulsar timing to fit for the position of the pulsar, as described in Ray et 
al. (2011). We measured 55 times of arrival (TOAs) based on 22-day integrations spanning 
the data set described above. The typical uncertainty on each TOA measurement was 4.7 
ms. Using Tempo2 (Hobbs et al. 2006), we fit the TOAs to a timing model including /, /, 
and /. With only these terms in the model, we observe very large residuals and the of 
the fit is very poor. This poor model fit means that the statistical errors in the fitted right 
ascension and declination reported by Tempo2 are unreliable. To get an estimate of the 
systematic error in the position fit, we use the following procedure. We added 5 so-called 
'WAVE" terms to the timing model to account for the timing noise using harmonically- 
related sinusoids (Hobbs et al. 2004). We then perform a fine scan over a positional grid 
around the location of source #20. Holding the position fixed at each grid point, we fit 
for the spin and WAVE parameters. The grid position with the lowest resulting for the 
fit is R.A. 20''21™29^683, Decl. +40°26'54.61" (J2000). This new timing position is 10" 
away from the one reported in Ray et al. (2011) which was based on 14 fewer months of 
data. Based on the map over the grid, we estimate the 95% confidence region of the new 
timing position to be an ellipse of dimensions 26"x 10", as shown in Figure |H The separation 
between the position of source #20 and the refined timing analysis position obtained here is 
14.7" and source #20 lies outside the 95% confidence region. 

A precise evaluation of systematic timing errors is complicated by the (erratic) timing 
behavior of the pulsar itself and the relatively short data span with respect to the 1-year 
modulation that is introduced by an incorrect position. That is, the position can be per- 
turbed by any component of the timing noise that appears to be a 1-year sinusoid. The 
magnitude of this effect is difficult to estimate because we have just one realization of the 
stochastic timing noise process in our data. To see the potential contribution, we have 
plotted an estimate of the timing noise contributions in Figure O Here, we assume that the 
parameters / and / are dominated by the secular spindown of the pulsar, while any higher 
order frequency derivatives and all of the WAVE parameters are dominated by timing noise. 
Clearly, this is a red-noise stochastic process. Note that a position error of 10" will introduce 
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a sinusoidal term with a 1-year period and an amplitude of 24 ms. To get another estimate of 
the systematic error, we fit a timing model including the pulsar position to three overlapping 
segments of data, the first half of our data span, the middle half, and the last half. These 
three fitted positions are separated by 7-9", giving another estimate of the systematic error 
resulting from timing noise. 

Based on these considerations, we adopt 10" as an estimate of the systematic error in 
the timing position. Combining, in quadrature, this systematic error estimate with the error 
estimated in Ray et al. (2011) results in the smaller, 10.3" radius, circular 95% confidence 
region depicted in Figure HI Source #20 lies within this region. 

In summary, when using pulsar timing to derive a position at the few arcsec level it is 
important to allow for low frequency (year timescale) timing noise. We have done this two 
ways: First, by adding WAVE terms to the solution and allowing these terms to vary when 
deriving an error in our grid search. Second, by time slicing the data and looking at how 
the derived position changes. As seen in Figure SJ each of these methods gave similar sized 
error regions that overlap. 

4. Search for an Optical Counterpart to the X-Ray Source 

As reported in Weisskopf et al (2006) and repeated in §5.3[ there are no cataloged optical 
counterparts for source #20. This is not surprising as, with the exception of the m ~ 16 
Crab, most optically detected pulsars have magnitudes > 25. In addition the field of PSR 
J2021+4026 is crowded and optical observations are further hampered by the presence of 
the 8^^ magnitude star BD+39 4152 (=V405 Cyg) one arcmin away to the east. 

We present here observations of the field taken on October 31, 2008 with the OPTIC 
orthogonal frame-transfer camera on the Kit Peak National Observatory, 3.6m, Wisconsin, 
Yale, Indiana, & NOAO (hence WIYN) telescope. The OPTIC camera, with a 10' field and 
plate scale 0.141"/pixel, allows improved image quality through "Orthogonal Transfer" (OT) 
rapid electronic guiding following motions of a reference star (Tonry et al. 2004). We used 
BD-l-39 4152 itself as the guide star and were able to correct at 50 Hz, collecting 3 x 180 s 
dithered exposures in r' and i'. These frames were subject to standard calibrations, except 
for the fiat field frames which were assembled by applying image shifts matching those of 
the OT guiding during the individual science exposures. The resulting image stacks have 
final PSF widths of 0.87" (r') and 0.62" {i') near the guide star; the PSF width increases by 
~ 30% toward the edge of the frame. We estimate that the frame is aligned to the Chandra 
coordinates with < 0.2" precision. 
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Figure |6] shows a portion of the OPTIC frames in r' and i' centered near the position 
of source #20. Note the secondary reflections of BD+39 4152 and the strong scattering 
background, especially in i' . Magnitudes were corrected to the SDSS photometry scale using 
observations of the calibration star Ru 149F (Smith et al. 2002). We measured the fluxes 
of the faintest detectable stellar sources in the vicinity of our target position and used these 
to estimate upper limits (~ 95% confidence) on the undetected optical flux for source #20 
of i' > 23.0 mag (the sensitivity is severely limited by scattered flux) and r' > 25.2 mag. 

Some of the diffuse emission toward the right (west) of the r' image is part of larger scale 
filamentary structure visible over several arcminutes. This is likely Ha/ [Nil] associated with 
the 7 Cygni SNR itself. We note that this remnant has been poorly studied in the optical. 
Mavromatakis (2003) described extensive diffuse line emission over a ~ 1° region, but found 
little filamentary emission and was not able to detect the very faint filaments seen in our 
data. No corresponding X-ray structure is seen, supporting the claim in Mavromatakis 
(2003) that the 7 Cygni remnant is dominated by low velocity shocks. Deep narrow band 
imaging to trace this structure could be useful in testing the connection, if any, between 
PSR J2021+4026 and the 7 Cygni SNR. 



5. The Other 43 X-Ray Sources in the Field 

5.1. Spectral Analysis 

There is a drop in the Chandra energy response above the mirror coating's iridium-M 
edges (~2 keV), thus any source with a substantial fraction of its detected photons above 
2 keV is indicative of a very hard spectrum. Figure [7] shows the X-ray color-color diagram 
for the 21 sources that have more than 15 source counts. The diagram comprises 3 bands: S 
(soft) covering 0.5 to 1.0 keV; M (medium) covering 1.0 to 2.0 and H (hard) 2 to 8 keV with 
T (total) simply the sum of S, M and H. The color ratios that comprise Figure [7] are given 
in columns 8 & 9 of Table [1] The X coordinate in Figure [7] (H-S) /T measures how hard 
the spectrum is and the y coordinate (M/T) measures how centrally peaked the spectrum 
is. Positive source counts requires data points to be inside the triangle, but background 
subtraction causes a few to appear slightly outside. Sources with (H-S)/T greater than 
0.5 are spectrally very hard and are likely background AGNs shining through the galactic 
plane. Only one of this group of sources, #42, has a cataloged optical counterpart (Table |5]). 
Sources with negative values of (H-S) /T likely have thermal spectra and are plausibly lightly 
absorbed foreground stars. Note that source #20 has the highest fraction of counts in the 
1-2 keV band. Two X-ray sources, other than source #20, have sufficient counts to warrant 
a spectral analysis: #3, the brightest in the field; and #32. 
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5.1.1. The X-Ray Spectrum of Source ^3 

We first fit the source #3 data using absorbed powerlaw, bbodyrad, and mekal mod- 
els. None of tliese models provided acceptable fits to the data, being 48.3, 57.2, 70.8, 
respectively, for 32 degrees of freedom. The only two-component+ model that provided an 
acceptable fit was a two-temperature mekal model with of 23.4 for 30 degrees of free- 
dom, further indicating that this is a foreground star. The results of our spectral analysis, 
following the procedures discussed in §2.2. H are in Table HI 

5.1.2. The X-Ray Spectrum of Source #32 

The data for source #32 are well-fit by an absorbed powerlaw model (x^ of 15.7 for 
19 degrees of freedom), but, as with source #20, the powerlaw index is very steep being 
3.4. In this case, however, neither a single-temperature mekal model nor a bbodyrad model 
provide as compelling fits (x^ of 46.9 & 23.9 respectively), although the latter is statistically 
acceptable. There is simply too much uncertainty to firmly classify this source on the basis 
of the X-ray spectroscopy alone. 

5.2. Temporal Variability 

The general paucity of counts also precludes a sensitive time-variability analysis for 
almost all these sources. Nonetheless, one of the three X-ray-brightest sources, #32, shows 
evidence for a significant temporal variation, suggestive of stellar coronal emission. The 
existence of both a likely 2MASS and USNO candidate counterpart (Tabled]) reinforces this 
interpretation in which case simple spectral models might not be expected to fit the data, 
as we have seen ( §5.1.2p 

5.3. Candidate Catalog Optical and Near-Infrared Counterparts 

We searched for candidate optical counterparts to the detected X-ray sources. We used 
HEASARC's BROWSE0 feature to search for cataloged objects within the 99%-confidence 
radius (egg) of X-ray source positions in Tabled! Table [5] tabulates results of a cross correla- 
tion of the X-ray positions of the Ciiandra-detected sources (column 1) with optical sources 



See 



http://heasarc.gsfc.nasa.gov/db-perl/W3Browse/w3browse.pl 
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(columns 2-7) in the USNO-Bl.O catalog (Monet et al. 2003) and with near-infrared sources 
(columns 8-12) in the 2MASS catalog (Skrutskie et al. 2006). 

5.3.1. USNO-Bl.O 

For fourteen (14) X-ray sources, we found a USNO-Bl (optical) source within the 99%- 
confidence radius egg of the Chandra position (Tabled]). Table [5] columns 2-4 list, respec- 
tively, the USNO-Bl right ascension RA, declination Dec, and RMS positional error cTq in 
the form (cXo (RA), do (Dec)). Column 5 gives the angular separation 6ox between optical 
and X-ray positions; column 6, the I-band magnitude. Column 7 estimates the probability 
Po{Sox, I) for a chance coincidence within the observed separation of an object as bright or 
brighter than the I magnitude of the optical candidate. We determined this probability 
from the I-magnitude distribution of the 893 USNO sources within 6' (slightly larger than 
the 8x8' Chandra field of view) of the X-ray pointing direction. We designate a potential 
optical counterpart to an X-ray source as a 'strong candidate' only if the sample impurity — 
i.e., probability of chance coincidence — Poi^ox,^) < 1%- AH the candidate USNO-Bl sources 
satisfy this criterion. 

5.3.2. 2MASS 

For nineteen (19) X-ray sources, we found a 2MASS (near-infrared) source within the 
99%-confidence radius egg of the Chandra position (Table [T]). Table E] columns 8 and 9 
list, respectively, the 2MASS right ascension RA and declination Dec, each with an RMS 
positional error cTj ^ 0'.'080. Column 10 gives the angular separation 6ix between near- 
infrared and X-ray positions; column 11, the K^-band magnitude. Column 12 estimates the 
probability Pi{6ix,Ks) for a chance coincidence within the observed separation of an object 
as bright or brighter than the Kg magnitude of the infrared candidate. We determined 
this probability from the i^^-magnitude distribution of the 1188 2MASS sources within 6' 
of the X-ray pointing direction. We designate a potential near-infrared counterpart to an 
X-ray source as a 'strong candidate' only if the sample impurity — i.e., probability of chance 
coincidence — Pi{6ix,Ks) < 1%. Eighteen (18) sources satisfy this criterion, the exception 
being source #44. Note that the 2MASS set of 18 strong candidate counterparts includes 
13 of the 14 USNO-Bl set of strong candidates f^ lSXTD . 

Table E] tabulates the 2MASS near- infrared photometry (columns 2-7) of the 18 strong- 
candidate optical (visible-near-infrared) counterparts to Ciaiidra-detected X-ray sources. 
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Examination of the near-infrared color-color diagram for all 2MASS sources within 6' of the 
pointing direction indicates those that are strong-candidate counterparts to X-ray sources 
are distributed as the field sources — i.e., as reddened main-sequence stars. Although most 
Galactic-plane 2MASS objects are normal stars, the objects identified with the X-ray sources 
need not be normal stars: For example, the X-ray emission may originate in an accreting 
compact companion. Figure [8] shows the X-ray hardness ratio of the 9 brighter X-ray sources 
versus the infrared color J-H of the corresponding strong 2MASS counterparts. The x-ray- 
soft and bluer infrared sources in the lower left corner of the figure are likely foreground 
stars. The x-ray-hard and reddened sources in the upper right hand corner may be X-ray 
binaries and/or AGN. 

5.3.3. WIYN observations 

The images we describe in Section [Hallow us to measure or limit the optical magnitudes 
of several other X-rays sources. A few sources were not covered in all sub-frames of the image 
stack and so were measured from individual exposures. Three additional optical counterparts 
to the X-ray sources were detected in the low exposure guide sector allowing improved frame 
registration. Table [7] gives the detected magnitudes and upper limits. 

6. Further Discussion and Summary 

Using the Chandra X-ray Observatory, we continued our search (Becker et al. 2004, 
Weisskopf et al. 2006) for possible X-ray counterparts to the intriguing 7-ray source now 
known as PSR J2021-I-4026. We found 44 X-ray sources in a field centered on the PSR 
J202H-4026 position, located along the line-of-sight toward the 7-Cyg SNR. Only one of these 
sources, #20, can reasonably and with high confidence be taken as the X-ray counterpart to 
the 7-ray source. 

There are a number of reasons supporting this conclusion. First and foremost, our X-ray 
source #20, is only 14.7" distant from the best-fitting 7-ray timing position. In addition, 
this separation is within the combined statistical and systematic errors on that position. 
There are also no other X-ray sources detected within 66" of the 7-ray source position to a 
Chandra source-detection limit of ~10~^^ erg cm~^ s~^ in the 0.5—8.0 keV bandpass making 
it highly unlikely that any of the other X-ray sources in the field are candidate counterparts. 
Furthermore, the spectrum of source #20 has a shape consistent with soft (log(Too) ~ 6.0 
to 6.5) thermal emission as expected from a young neutron star though perhaps somewhat 
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higher than expected from the spin-down age of 77000 yr estimated for PSR J202 1+4026. 
There is also a hint of extended diffuse X-ray emission in the vicinity of source #20 that 
may be an associated PWN. 

With source #20 as the counterpart, we infer F-y/Fx ~ 1.1 x 10^, not atypical of young 
isolated neutron stars (e.g., Becker 2009). If source #20 is not the X-ray counterpart, the 
flux ratio is at least 30 x larger, which would be substantially larger than the observed ratio 
for other 7-ray pulsars. 

A similar argument using the optical data also supports source #20 as the counterpart: 
our r' > 25.2 limit implies a lower limit of Fx/ Fy ~ 250, with some uncertainty due to 
extinction. This is already larger than the maximum value for X-ray binaries (~15) or BL 
Lacs (~100) and is approaching typical values for isolated neutron stars (~ 10'^^'^) (e.g. 
Schwope et al. 1999). Thus, based on the X-ray/optical evidence alone, source #20 is likely 
an isolated neutron star and is the likely counterpart for PSR J2021+4026. 

Finally, the X-ray spectrum has a shape consistent with the soft thermal emission ex- 
pected from a young neutron star. This emission likely represents a fraction of the stellar 
surface heated by back-flowing particles generated by magnctospheric activity (eg. Harding 
& Muslimov 2001). At present the fitted parameters suggest that this thermal component 
implies a relatively large distance, ~ 6 kpc, incompatible with an association with SNR 
G78-I-2.1. However, the fits also indicate a complex spectrum with at least two components; 
much higher S/N data will be needed to extract strong spectral constraints. Of course, a 
heated polar cap suggests that sensitive observations should also be able to detect X-ray 
pulsations at the 265 ms spin period, the definitive test of the counterpart's association. 
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Table 1. Chandra X-ray sources on ACIS-S3 in the Fermi-LAT PSR J02 1+4026 field. 
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1 


.4 


26 


0, 


.34 


0, 


.69 


0.14 ± 


0, 


,12 


0.71 ± 0.09 


23 


20 


21 


31, 


.889 


40 24 25. 


,51 


1 


.4 


11 


0. 


.39 


0, 


.79 










24 


20 


21 


32 


.669 


40 28 21. 


,38 


2, 


.2 


11 


0. 


.50 


1. 


.01 










26 


20 


21 


32 


.906 


40 24 20. 


.88 


0, 


.8 


14 


0. 


.33 


0, 


.66 










26 


20 


21 


33, 


.031 


40 23 0, 


.62 


2, 


.3 


11 


0. 


.52 


1 


.06 










27 


20 


21 


33, 


.650 


40 29 8, 


.97 


2 


.8 


48 


0, 


,39 


0, 


.78 


0.80± 


0, 


11 


0.25 ± 0.09 


28 


20 


21 


34 


.097 


40 25 26 


50 


1 


.3 


12 





38 





.76 










29 


20 


21 


34, 


.559 


40 23 19, 


.16 


2, 


.0 


33 


0, 


,37 


0, 


.73 


0.86 ± 


0, 


,10 


0.09 ± 0.07 


30 


20 


21 


36 


.268 


40 28 35. 


.85 


2, 


.0 


16 


0. 


,43 


0, 


.87 










31 


20 


21 


36 


.486 


40 28 13. 


.57 


1 


.9 


6 


0. 


,55 


1 


.10 










32 


20 


21 


37, 


.679 


40 29 68. 


.09 


3 


.6 


248 


0. 


33 


0, 


.67 


-0.21 ± 


0. 


,05 


0.66 ± 0.03 


33 


20 


21 


38, 


.401 


40 29 35. 


.49 


2, 


.8 


53 


0. 


,38 


0, 


.76 


0.74 ± 


0. 


,11 


0.22 ± 0.08 


34 


20 


21 


38, 


.431 


40 24 42, 


77 


1 


.7 


94 


0, 


32 


0, 


.64 


-0.45 ± 


0. 


.07 


0.43 ± 0.05 


35 


20 


21 


38, 


.579 


40 24 14 


88 


1 


.3 


7 





43 


0, 


.86 










36 


20 


21 


39, 


.214 


40 27 9, 


.90 


1 


.0 


9 


0, 


,36 


0, 


.73 










37 


20 


21 


40, 


.083 


40 24 9, 


.43 


2, 


.6 


11 


0. 


,55 


1. 


.11 










38 


20 


21 


43, 


.107 


40 23 53, 


.76 


2, 


.6 


48 


0, 


,38 


0, 


.76 


-0.08 ± 


0, 


.10 


0.72 ± 0.07 


39 


20 


21 


44 


.548 


40 29 34 


34 


3 


.3 


19 


0, 


,55 


1 


.10 










40 


20 


21 


47, 


.294 


40 24 54 


84 


3 


.2 


28 


0, 


47 


0, 


.95 


0.99 ± 





16 


0.01 ± 0.12 


41 


20 


21 


47, 


.584 


40 26 57, 


50 


1 


.9 


18 


0, 


41 


0, 


.82 


-0.53± 


0, 


44 


0.20 ± 0.13 


42 


20 


21 


52, 


.529 


40 25 7, 


.72 


4, 


.2 


66 





44 


0, 


.87 


0.79 ± 


0, 


,09 


0.16 ± 0.07 


43 


20 


21 


57, 


.170 


40 26 24, 


.36 


7, 


.8 


29 


0, 


,93 


1 


.87 


0.94 ± 


0, 


,23 


0.15 ± 0.18 


44 


20 


21 


67, 


.761 


40 26 46. 


.78 


6, 


.8 


106 


0. 


,46 


0, 


.92 


-0.42 ± 


0. 


,12 


0.31 ± 0.04 



Table 2. Fits to single models with absorption. 



parameter 


powerlaw 


bbodyrad 


nsa" 


nsmax-1260° 


nsmax-130190'' 




l.Oo 


U. ( y 


1 flQ 

i.Uo 


l.Uo 


i.Uo 


C 


i n on 1 n oo^ 


n 1 1 n 1 q\ 
t^-U.lo,+U.le5j 


t^-U.iy,+U.zi ) 


('nin inoi\ 
(-U.iy,+U.zi ) 


/ n 1 Q in oi ^ 
(^-U.io,+U.zl } 


pl-index 


4.86 










^yl index 


(-0.48,+0.55) 










log(Toc, 




6.52'' 


6.16 


6.12 


6.10 






(-0.05,+0.04) 


(-0.08.+0.07) 


(-0.08,+0.08) 


(-0.08,+0.08) 


norm 


1.06 xlO-4 


0.80 


1.8x10-8 


0.56 


0.83 


^ norm 


(-0.32,+0.51)xl0-4 


(-0.37,+0.83) 


(-1.2,+4.5)xl0-9 


(-0.38.+1.44) 


(-0.19,4-0.21) 


D(kpc)'= 




174. 


23 


16 


13 






(-52,+63) 


{-4,+17) 


(-8,+12) 


(-1.4,+1.8) 


Rem (km)-'' 




0.1 


0.9 


1.2 


1.5 




23.3 


30.6 


29.8 


28.5 


27.7 


degrees of freedom 


22 


22 


22 


22 


22 


Flux xlO^'* (ergcm"-^ s~^) 


43 


4.9 


8.0 


8.9 


2.0 



Note. — 

" The mass, radius, and magnetic field were fixed at 1.358Mq, 12.996km, and 1.0 x 10^* gauss respectively. 

^ The mass and radius of the neutron star were chosen as for the nsa model so that the input, the gravitational redshift IH-z, 

was fixed at 1.15. 

Uncertainties based on considering only 1 interesting parameter, i.e. the bounds indicated by the minimum x^-|-l. 

Assumes M/R=l. 358/12. 996 Mg/km. 

Derived from the different normalizations assuming a NS mass and radius 1.358Mq, 12.996km. 
f Rough estimate of the size of the emitting region by scaling the distance to 1.5 kpc. 
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Tablc 3. Fits to dual models with absorption. 



parameter 


bbodyrad 


nsa" 


nsmax-1260'' 


nsmax-130190'' 




16.38 


16.24 


16.36 


16.38 


degrees of freedom 


20 


20 


20 


20 


f statistic 


4.23 


8.34 


7.03 


6.93 


probability 


0.029 


0.0023 


0.0049 


0.0052 


IlH 


0.76 


0.97 


0.98 


0.98 




(-0.16,+0.19) 


(-0.18, +0.22) 


(-0.19,+0.23) 


(-0.18,+0.23) 


log(Too 


6.41 


6.01 


5.98 


6.00 




(-0.08,+0.07) 


(-0.12.+0.11) 


(-0.13,+0.11) 


(-0.14,+0.10) 


D'*(kpc) 


91. 


6.8 


6.0 


5.7 


CD " 


(-15, +39) 


(-4.7,+9.5) 


(-4.3,+9.4) 


(-4.4,+6.2) 


Rer„(km)9 


0.2 


2.8 


3.3 


3.4 


r 


1.2 


1.10 


0.73 


0.73 




(-1.5,+1.2) 


(-1.6,+1.3) 


(-2.10,+1.44) 


(-2.2,+l,5) 


PL norm 


1.5 


1.26 


0.69 


0.68 


O'norm ^ 10^ ^ 


(-1.5,+6.0) 


(-1.13,+5.88) 


(-0.57,+1.05) 


(-0.56,+7.02) 


Flux-f xlO" (ergcm-2s-i) 


8.4 


15.8 


16.3 


16.8 



Note. — 

" The mass, radius, and magnetic field were fixed at 1.358M0, 12.996km, and 1.0 x lO^^ gauss 
respectively. 

^ The mass and radius of the neutron star were chosen as for the nsa model so that the gravitational 
redshift 1+z was fixed at 1.15. 

Derived from the different normalizations assuming, where necessary, a NS mass and radius of 
1.358M0, 12.996km respectively. 

Uncertainties based on considering only 1 interesting parameter, i.e. the bounds indicated by the 
value of the parameter in question at the minimum x^+l. 
f Unabsorbed flux. 

^ Rough estimate of the size of the emitting region by scaling the distance to 1.5 kpc. 



Table 4. Fit to a two-temperature mekal model for source #3. 



nn 


feTi 


kT2 


normi 


norm2 


(1022 cm-3) 


(keV) 


(keV) 


xlO* 


xlO^ 


0.44 


0.20 


1.03 


1.37 


2.36 


(-0.13, +0.11) 


(-0.02, +0.04) 


(±0.07) 


(-0.9. + 1.9) 


(-0.32, +0.25) 



Note. — 

Uncertainties based on considering only 1 interesting parameter, i.e. the bounds 
indicated by the value of the parameter in question at the minimum 



Table 5. Candidate cataloged counterparts to X-ray sources in the PSR J2021-I-4026 field. 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



USNO (optical) candidate counterpart 2MASS (infrared) candidate counterpart 

X-ray Source RA(J2000) Dec(J2000) ctJ 6ox I Po{6ox,l) RA(J2000)'' Dec(J2000)'' Si^ Pi{6ix,'^a) 

h m B o / // // // jjja^g % h m s o , „ „ jjjg^g % 



1 


20 


21 


09.036 


40 


27 


03.51 


(0.04,0.07) 


0.24 


18.58 


0.36 


20 


21 


09.035 


40 


27 


03.30 


0.12 


13.625 


0.18 


3 


20 


21 


11.433 


40 


28 


05.01 


(0.04,0.10) 


0.18 


13.15 


0.07 


20 


21 


11.424 


40 


28 


04.71 


0.18 


11.155 


0.01 


5 






















20 


21 


13.324 


40 


26 


02.91 


0.24 


14.501 


0.27 


6 






















20 


21 


13.596 


40 


25 


55.80 


0.24 


14.175 


0.21 


8 


20 


21 


14.343 


40 


25 


20.63 


(0.04,0.03) 


0.24 


14.91 


0.12 


20 


21 


14.348 


40 


25 


20.37 


0.12 


12.064 


0.04 


9 


20 


21 


06.967 


40 


25 


17.44 


(0.04,0.04) 


0.24 


13.71 


0.10 


20 


21 


16.967 


40 


25 


17.17 


0.06 


12.559 


0.06 


10 


20 


21 


19.557 


40 


25 


33.02 


(0.04,0.03) 


0.42 


18.19 


0.30 


20 


21 


19.549 


40 


25 


32.59 


0.36 


13.869 


0.20 


11 


20 


21 


20.702 


40 


24 


02.14 


(0.04,0.07) 


0.30 


14.58 


0.12 


20 


21 


20.700 


40 


24 


01.97 


0.12 


12.249 


0.05 


15 


20 


21 


25.749 


40 


27 


44.39 


(0.04,0.08) 


0.30 


17.21 


0.22 


20 


21 


25.767 


40 


27 


44.06 


0.18 


13.722 


0.18 


17 


20 


21 


28.657 


40 


24 


15.97 


(0.04,0.02) 


0.12 


19.07 


0.35 


20 


21 


28.679 


40 


24 


15.50 


0.42 


15.083 


0.44 


18 






















20 


21 


29.779 


40 


24 


55.53 


0.42 


15.002 


0.35 


22 


20 


21 


31.425 


40 


22 


56.71 


(0.04,0.16) 


0.54 


17.47 


0.18 


20 


21 


31.422 


40 


22 


56.41 


0.42 


13.188 


0.09 


32 


20 


21 


37.566 


40 


29 


58.32 


(0.04,0.07) 


0.30 


14.54 


0.09 


20 


21 


37.575 


40 


29 


57.90 


0.18 


13.302 


0.09 


34 


20 


21 


38.410 


40 


24 


42.97 


(0.04,0.08) 


0.30 


13.11 


0.08 


20 


21 


38.427 


40 


24 


42.63 


0.12 


11.476 


0.02 


37 






















20 


21 


40.056 


40 


24 


09.63 


0.36 


14.465 


0.65 


38 


20 


21 


43.117 


40 


23 


54.22 


(0.04,0.08) 


0.48 


13.68 


0.11 


20 


21 


43.116 


40 


23 


54.05 


0.30 


11.967 


0.04 


39 


20 


21 


44.515 


40 


29 


34.03 


(0.04,0.03) 


0.48 


18.00 


0.55 


20 


21 


44.538 


40 


29 


33.44 


0.90 


14.003 


0.45 


41 


20 


21 


47.628 


40 


26 


57.94 


(0.04,0.32) 


0.66 


17.33 


0.24 




















42 






















20 


21 


52.555 


40 


25 


08.16 


0.54 


13.863 


0.25 


43 






















20 


21 


57.240 


40 


26 


23.10 


1.50 


14.723 


2.28 



Note. — " USNO RMS positional uncertainty in each axis (RA, Dec) 
2MASS RMS positional uncertainty = 0'.'08 per axis 
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Table 6. 2MASS infrared photometry of strong-candidate counterparts to X-ray sources. 

















1,4) 


X-ray Source 




J 






J-H 




H-Ks 


1 


15 


.583 ± 


.043 


1 


.436 ± 0, 


,068 


0.522 ± 0.073 


3 


11 


.927 ± 


.018 





.675 ± 0, 


,028 


0.097 ± 0.027 


5 


16, 


.031 ± 0, 


.068 


1 


.285 ±0, 


.107 


0.245 ± 0.133 


6 


15, 


.772 ± 0, 


.056 


1, 


.195 ±0, 


.087 


0.402 ± 0.099 


8 


13, 


.113 ±0. 


.021 


0, 


.698 ± 0. 


.030 


0.351 ± 0.030 


9 


13 


.045 ± 0, 


.021 





.387 ± 0. 


.030 


0.099 ± 0.032 


10 


15, 


.278 ± 0, 


.040 


1, 


.177 ±0. 


.059 


0.232 ± 0.073 


11 


13 


.154 ± 0, 


.022 


0, 


.713 ±0, 


.033 


0.192 ± 0.034 


15 


15, 


.133 ± 0, 


.038 


0, 


.991 ±0, 


.064 


0.420 ± 0.074 


17 


16, 


.021 ± 0, 


.089 


0, 


.945 ± 0, 


.121 


-0.007 ±0.190 


18 


16, 


.046 it 0, 


.136 





.442 ± 0. 


.164 


0.602 ±0.180 


22 


14 


.761 ± 0, 


.027 


1, 


.238 ± 0. 


.041 


0.335 ± 0.047 


32 


13, 


.887 ± 0, 


.021 


0, 


.534 ± 0. 


.031 


0.051 ± 0.043 


34 


12, 


.045 ± 


.021 





.517 ±0, 


.030 


0.052 ± 0.028 


37 


15, 


.822 ± 0, 


.064 


1, 


.119 ±0, 


.093 


0.238 ± 0.124 


38 


12, 


.582 ± 0, 


.021 


0, 


.480 ±0, 


.030 


0.135 ± 0.028 


39 


15, 


.477 ±0, 


.108 


0, 


.974 ± 0, 


.129 


0.500 ± 0.123 


42 


15, 


.995 ± 0, 


.062 


1, 


.307 ± 0. 


.097 


0.825 ±0.101 


43 


16 


.249 ± 0, 


.075 


1, 


.280 ± 0. 


.128 


0.246 ±0.160 



Table 7. WIYN Optical MagnitTidcs 



Source 


i' 


r' 


Source 


i' 


r' 


Source 


i' 


r' 


Source 


i' 


r' 


1 


18.89 


21.04 


2 


~24.4 


>" 


3 


15.52 


15.77 


4 


17.96 


19.95 


5 


19.30 


21.36 


6 


19.02 


19.02 


7 


> 


> 


8 


15.63 


17.07 


9 


15.15 


15.50 


10 


18.39 


20.39 


11 


15.63 


15.75 


12 


> 


> 


13 


22.40 


> 


14 


> 


> 


15 


17.94 


19.82 


16 


> 


> 


17 


19.37 


21.83 


18 


18.13 


20.15 


19 






21 


> 


> 


22 


17.82 


19.94 


23 


20.92 


23.29 


24 


~ 24.4 


~ 25.5 


25 


> 


> 


26 


> 


> 


27 


22.45 




28 






29 


> 


> 



^Magnitude limits: > = > 24.2 (i'), = > 25.3 (r') 
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3500 4000 4500 



Fig. 1. — The field showing the most recent Chandra observation, ObsID 11235. For this 
figure a nearest neighbor smooth has been apphed. The rotated square shows the extent of 
the ACIS S3 chip, i.e., the region searched. The color bar on the right shows the number of 
(smoothed) counts detected in a pixel during the observation. The numbers on the left and 
bottom show the ACIS coordinates in pixels. Sources are numbered in order of increasing 
right ascension. See Table 1 for source X-ray properties. 
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H — h-^■ 



Energy (keV) 

Fig. 2. — The spectrum of source ^20 with background subtracted and fit to a nsa+powerlaw 
model. The dotted hne in the upper panel is the nsa component, the dashed line is the 
powerlaw component and the solid line is the two components combined. The lower panel 
shows the contributions to x^. 
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Fig. 3. — Significance of pulsation detection, as measured by tlie test, versus source 
number. Tfie data were binned into 20 pulse phase bins and fit to a constant fiux model 
(no pulsations). The worse the fit (higher x^), the more likely the X-ray source is the 7-ray 
pulsar counterpart. 
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Fig. 4. — Same as Fig. [T]but now just showing the region around source #20. The large 
circle denotes the most recent LAT imaging position 95% confidence error circle (2FGL, 
Abdo et al. 2011). The small ellipse is the 95% confidence timing ellipse from the current 
work using WAVE terms to estimate the impact of timing noise. The small circle is the 
95% confidence region obtained by combining the timing solution determined by Ray et al. 
(2011) added in quadrature with our new (10") estimate of the systematic error. See text 
for further details. 
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l_i I I I I I I I I I I 

5.48x10^ 5.5x10^ 5.52x10^ 5.54x10^ 5.56x10^ 

Time (MJD) 

Fig. 5. — Illustration of the magnitude of the timing noise observed in PSR J2021+4026. The 
dotted (blue) curve shows the timing residuals that can be attributed solely to timing noise. 
This is the sum of two components: The thin upper (red) curve is the contribution from 
the polynomial terms / and / , neither of which can be attributed to the secular spin down 
of the pulsar. The lower (green) curve is the combined sinusoidal WAVE components (see 
text). Note that the pulsar period is 0.265 seconds, so timing noise causes several complete 
phase wraps over this interval. 
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Fig. 6. — Portions of the WIYN 3.6m/OPTIC r' (left) and i' (right) images of the field of 
PSR J2021+4026. Reflected light and scattering from the magnitude star BD+39 4152, 
used here as a guide star and located to the left of each image, dominates the fields. The 
position of Chandra source #20 is marked with a circle of radius 2 arcsec. There is no 
detection visible in either band with upper limits as given in the text. 
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-1 -0.5 0.5 1 

(H-S)/T 

Fig. 7. — The X-ray color-color diagram for the 21 X-ray sources that have more than 15 
source counts. The Chandra band is divided into 3 bands with S (soft) covering 0.5 to 
1.0 keV, M (medium) covering 1.0 to 2.0 and H (hard) 2 to 8 keV. The Total, T, is the 
sum of S, M and H. The solid curves within the triangle represent powerlaw spectra with 
photon indices of 1, 2, 3 and 4. The dotted curves correspond to hnes of constant absorbing 
columns of {bottom to top) 1 x lO^o, 1 x lO^^, 2 x 10^1,5 x lO^^, and 1 x lO^^, respectively. 
The triangle encloses the physically-meaningful range of colors. See text for the discussion 
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Fig. 8. — X-Ray hardness ratio (H-S)/T defined in the caption to figure [7] versus the near 
infrared color J-H for X-ray sources with 15 or more counts and hkely 2MASS candidate 
counterparts. 
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